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ABSTRACT 


It is the purpose of this thesis to calculate the date re- 
quired to construct enthalpy-—entropy diagrams for propéilant 
Systems and to investigate their usefulness when applied to the 
cetermination of rocket motor performance parameters, 

Tiagrams which cover 4 range of mixture ratics bracketing 


the stoichiometric value are included for the fcllowing systems; 


(1) Liquid Ammonia and RFNA 
(2) Liquid Ammonia and Liquid Ozone 
(3) Liquid Hydrogen and Liquid Ozone 


Graphical perfcrmance parameter values obtained for the 
Jt diy tiwnthest mixture of liquid amnonta anc HFA differed less 
then two tenths of one percent from the analytical values cal- 
culated by the Jet Propulsion Laboratory, 

Performance values ere determined for chamber pressures of 
GOOG PSIA and 300 PSIA and for exhaust pressures of 14,7 P&IA, 
7.05 PSIA, and .147 PSIA. 

Because of the extensive calculations required to construct 
an h-s diagram, the graphical method of analysis is net eco- 
nomical unless a comprehensive study of « propellant is desirea, 
However, it performance velues are required over a wice range of 
boundary conditions, the h-s diagram .permits the presentation 
of a lerge amount cf data in 4& concise ena readily useable 


Lorm. 
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INTOGDISTION 


Cne of tre chief cisudvanteges of Cetermining rocket motcr 
performance parameters by the analytical methods based on the 
use of isentropic flow e,uctions is that one set of calculate 
icons yields results which are applicable for only @ very limited 
range of physical boundcery conditions, For exauple, values céel- 
culated for the expunsicn to an ambient exhaust pressure of one 
atmosphere cunnot be applied to the evaluution of u moter which 
is to propel the last stige of & multiple-step rcecket, Als, 
twe sets of calculations sre required te obtain the performance 
fcr both equilibrium anc ccunstaunt composition fiow, 

One sclution to these problems might be the graphical pre- 
sentiticn of the fundamental auta requirec to calculate the 
thecretic.l performence parameters for rocket motors in Such é 
vey aS to be useful over a side renge of operating conditions. 
This technique of using plotted deta has been wicecly used in 
other ficslais #hén aniulytic.l methods have preven te be cumber- 
some, A notable excmple is the grephical analysis of high speed 
agirpline performance, | 

it is the purpese ot this report to prcevice enthalpy- 
entropy Giagrams for the procucts of combustion of several pro- 
prllent systems, and to investigate their usefulness when 
appliei te the determination of tne theoretical. performance of 
reeket motors, 

wiweac infcrme ticn rewuired rer the constructicn of the 
h-s diugrums is bye Gisektcd) Compan Cle of tne reacticn pro- 


ducts at avpropriate temperatures and pre: sures, The numuber of 


thes§’ component @ahcula@tions olticn must Be aad tc comstruct . 
given diavrem is greater then the nunber revuire.s fcr tre sane- 
lytical calculation of the vericreumes peramoters fcr cone set 
cf bouncary conditions, but the diagrim parmits & much more 
corprehensive sbucy ef the spesitieca propelient system. Also, 
onee the basic slugram is complatec, tne effort require: to 
plot & me orefékure lings, cac thus peraat the Gvaluation at 
tie system [Tr an entirely vuilferent oxhvast contitien, le net 
erat, 


intnolpy—entropy cdagrems for a rense ol mixture ratic 


which brackets the stcichicuetric value are constructed Tur the 


a 


Maan NI, wey » eer ee sere Ls opellent systems. The Ai wanks rv dS 


system was selecte? in orgor th.t narameter velues obtained 


grephicallv coule be compare. with the enmulyticsl] véines teint 


a! 


eéleulutea for this systen by the Jet Prcpulsicn buboratay. - 
The Nit. -O » ei bes Tp On, cy8tems were chcsen for two reasons: .irst, 


because of the current interest in Lli.wid o#ume «8S gn G.izent; 


ane seeomd, because thete syetems represent, retp@etively, ar 


C 


intervaciate anc a wery hct system, while the Ni -nlis syster. 
ws 2 relmively ctl cmé@, A> 4 1@hult, CHG camplese bee ci 


cn 


ct“ 


. 7 i. : wv 
baygreis covers & tewp.ratur® renge wp 6 S220 #.. Tie ree 
for the int@rest in liwuis dime if tne lurge® negative Wel 


ht tS Bet ct tormetion. 
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Vives into the rollo@ing pirtss 


(1} Liseucsien of assumption. ani. Introuncticn 


Uti! 
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Trometers. 


aout 
cf 


» (4) Sample Culeuletion. 
(3) Diseussion of Hesuits. 


(4) Conclusions ang Kecommenduticns. 


ig gare (1) the assumptions moace in caéleuleting the thermo- 
Gynemic date useu in the construction of the hes Gingrums are 
emphasized, and an attempt is meade te cleurly cefine the ref- 
erence levels for all values used. The as-umpticns which rust 
be considered when the diagrams ere used to ealculate perfor- 
mance are «also included in croer that the Limitetions cf the 
graphical method be understood, 


x 


The sémple calculuticn in pert (7: ig cerried cut in 
sufficient detail that the reacer should be eble te readily 
check any desired velue by following the exact procedure used 
in the criginul computation. Because of the wide range of ten- 
peratures und pressures fcr which minor ccmponent calculations 
were made, certain trencs in the computed velugs were noticed 
which probebly ¥ouLld not show up in a single anulytical cal- 
enulation, Metheds for utilizing these trends when mexing en ex- 
tensive set of computeticns are vkpliuined, 

It should be noted in cennecticn with parts (3) and (4) 


that the diseussion of resuits and the conclusicns «crasn from 


t 


this snelysis are based upen the -tucay of theoreticel maximum 
values of the performince parameters, enc thet no vu» perimental 
data Wes available for comperison, For this reesen, the ceon- 
struction of the cisgrauns themlelves is euphasizeé rather than 


the particular values .erivec frum then. 
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Number «ff moLer of E. 


+ 


in the recetion procucts 
Number of moles of HyG in the rusction products 
Tunbiur of moles cf CC in the reaction products 
Wunber of moles of Cl.. in the reacticn products 
wuniber of moles cf Ge in the rexetion products 
Rumber of moles of Ne in the reeetion prouucts 
Number of moles of HO in the reacticn procucts 
umber ef moles of CH in the renetion products 
Numoer of moles of EH in thé réelLotitm produdts 
Number of moles of QO in the reacthiem products 
Number cf gram items of carbon in the rcactznts 
Number of grem «tome cf hydrogen in the reactants 
Nunber oF grem stoms of nitrogen in the resetents 


Number of gram utumse GL oaygen in tae reccbtents 


“quilibrium constent busca cn the ratios of partisl 
eguiliprium constant busec om the ratios cl nolbe 


Particulsr equillbrium cunztants based cn the ratios 
ef mole treeticns of the component geses im given 
reucticns 

Total number cf acles cf the reacticn breducts at TK 

Huxber of moles of the ip, eccmpenent of the reacticn 


procucts at TOK 


Pe 


Pe 
TOK 


BA; 4 20K 


Teuperature of reectiin prouucts in the rocket sctor 
chamber before expansicn through the nozzle 

Temperature of reecticn products after expuensicn 
through the noz:le 

Chamber pressure in peuncs per syucre inch ebsolute 

fxhcust pressure in puunds per squere inch «bsolute 

Sensible heat of the 14, ccmpcnent cof the reaction 

products at TOK as referred to the cefined 
reference stcte 

Total enthelpy cf the rezction preducts at TOX whose 
reference state is the e.uilibrium composition at 
5L0CK 

Heat relbesed in keals by a combustien resecticn at T°K 

He.t of formation cf the ig, component of the reuction 
prouucts at SUC 

Number of grams of propellent mixture 

Totel enthelpy per gram of the rezecticn procucts at 
TOK whose reference state is the eguilibrium cca- 
position at 600% 

Total absolute entrcpy of the reaction products at 
T°K and a pressure of p atnuospheres 

Absolute entropy of the dey ecmponent ef the :eection 
products at T°K and 2 pressure cf une atmosphere 

Absolute entropy per gram of the reietion prouucts 


at TOK and a pressure of p etmespheres 


ry 


RENA 


Universel gas ccnotant, equal to 1.266 czlorics per 
grem mole per uegree Kelvin 

iufiective exhaust velocity in feet per second 

mechanical e@.uivelent cf heat, equal to <«.186 tol 
ergs per keal 

speeific im ulse in peuncs of thrust per pound of 
propellent burned per seecnd 

Acceleration of gravity ( arbitrarily fixed et 3.16 
feet per second in this report) 

vixture ratio. The ratio of the weight of cxidinmt 
to the weight of fuel in & propellent systen 


Red fuming nitrie ecid (nitric acid with 6.&* of NUo ). 
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PART I: Dit~CUSsIGN CF ...-cUMPTIUNS iN LNTHOLUCTIUN CF 


The assumpticns upon uhich this enelysis is based aust 
ineluce both those which «ere involved in the ecnstructicn of 
the diagrans themselves, and uny aaciticnal ones thet are re- 
“juired when the diagrams ere used in the culculation cf rocket 
motor performance parameters, Tne first group includes the 


following: 


(1) The ecuation of state 1s pve RT. 

Note: The source of the enthslpy and entropy tables 
usec is the report on Aap eR FF end the counct- 
icn cf State of the Propelient Gases" by Hirshfelder, 
NeClure, and Curtiss, reference (2). The basic equat- 
ion of state used in this report is a van der faals 
e.uaticn for which correcticns have been made for the 
overlapping of here spheres, However, this corrected 
ecuaticn Wes not usec in the caleulaticn of enthelpy 
and entropy velués, ana the perfect gas lew must be 


assumed, 


(2) There is no Jcule-Thompson effect. That is, for a 
constunt temperature, enthalpy does not vary as the pressure 


changes, 
(5) There is no enthulpy change on mixing. 


(4) The reference stete for all enthalpy values plottea . 


F ant 


Le the ecvilibrium c@Ggcittige ca FAC > wee prescure of cone 
etacsphere cf the gis ai: ture contzining the reacticn froaucts, 
it which thernocrnamic stite the enthalpy of the mixture is 
consicered to be zero. It sheulé be noted thet fcr all practi- 
cl purposes this is eyguivelent tc saying thet at the reference 
state mincr components have disiuppeured, 

This qualification upon the reference state that there be 
no mincr components present is not readily apparent because the 
tabulated values of enthslpy are all zero at ouo°K, However, if 
we consider the e-uilibriun expansicn through a nozzie of a 
specified prcepellent mixture to a tenpereture at which ne sig- 
nificant amcunts of the miner components remein, it is cbvivuus 
that the heat which has uppecrec aue te the recor.binetion of 
these minor ccempoments represents an inerease in energy which 
1s availuble te precuce un increate in the velecity cf the pro- 
pellent mixture. Thus, the enthilpy value of 4 point on the | 
hes ulugram ust represent net only the -ensible hest cbove 
SUC°K eonteined in the yas misture at the temperature corres — 
ponding to the point, but slso tue heat which is re ,ulrea to 
shift the equilibrium from that at the reference state tc that 
at the thermodynamic state representeca by the point. This 
Latter increment of heat is ercetly that which in an. equili- 
prium expansicn is muce avelisble for increased performence, 

In & constent ccnpesiticn expensicn the entire enthalpy 
Grop is cue tc the chenge in the sensible heat, and in proba- 


bly no actual expension is e.ullibrium ever ccapletely 


-. a . 


maintsjinec. In ganeral, the enth.loy change for an S24 puineder: 
« jf ; 

will be somewhat srecter than the change in sensible hect .ue 

te the Temoereture drop but less than the sum of the sensible 


heet chenge and the heat which would be rejected from the Au 


©} 


mixture if eauilibrium were maintained, 

It is necessary to define some reference level for enthal- 
py values because there is no thermodynumie state at which we 
een say that the enthalpy of all substances is zero, For this 
rezson, it is not possible to determine ubsolute enthalpy 
values, and any value used must be referred to some reference 


state. 


(5) The tabulated entropy velues used are absolute values 
of entropy of the gas components calculeted ct a temperature of 
T°K and a pressure cf cne atmosphere and considering the com- 
ponents to be pure perfect gzses. 

Because of the Third Law of Thermodynamics, which states 
that the entropy of most crystals at the sbsclute zero may be 
taken as zero, it is pessible to define and to Getermine ab- 
sclute entropics. Therefore, it is not nsecessury to aefine a 
reference state for entropy values in the same sense thet one 


was defined for enthelpy values. 


(6) All miner cumponents ure considered except atomic 
nitrogen (N). None cf the adiabatic flame temperatures for the 


HWHs-RFNA and NH.-O; systems greatly exceeas SuICTR, and fer 


- 3D w= 


this temperiture renge previous calculations have shown this 


uScupption to be justified. 


(7) Over a temperature range ef 100°K the variation of 
temperature dependent quentities such as gas mixture conposit- 


ion with temperature is assumed to be Linear. 


(8) The equipertition of energy among the electronic, 
vibrational, and rotetional energy levels is assumed to be in- 
stantaneous during the combustion and expansion processes for 


both equilibrium and constant composition flow. 


(9) The pressure lines on the h-s diagrams which repre- 
sent equilibrium flow sre calculated assuming that the compos- 
ition of the products of the combustion reaction changes con- 
tinuously to maintain tre equilibrium composition corresponding 
to the prevailing thermodynamic state, and that the energy re- 
leased by this equilibrium shift is completely available for 


conversion into energy of motion, 


(10) The pressure lines on the h-s diagrams which repre- 
“sent cconstint composition flow ere calculated essuming that the 
composition of the products of the combustion reaction remains 


constent end equal to the equilibrium composition at the adia- 


betic flame temperature T,. 


In using the h-s diagrans to calculate rocket motor per- 
formence parameters it is necessary to make the following 


adcitional assumptions: 


eee |) ee 


(11) The rexeticn of the propellent mixture proceeds te 
completion in the rocket motor chamber, and sufficient time is 
availeble for the equilibrium concentration of all components 
to be estublished at the adlebatic flane temperature determined 


by the mass and heat balsnce, 


(1E) Viscous effects are neglected. That is, all flew is 


ecnsidvred te be isentropic. 
(13) The gas flow through the nozzle is steady, 
(14) The gas flow leaving the nozzle is axial. 


(15) The velocity.of the gas mixture in the combustion 
chamber is assumed to be negligible in compérison with the 


effective exhaust velocity. 
(16) Gravity effects sere neglected, 


(17) The nozzle is always correetly expanded for the 
ambient exhaust pressure. If this pressure varies, the éssump- 


tion of & "rubber nozzle" must be mace, 


> 
7 Ros _ 


The performence parameters which are considered in this re- 
port are: | 

(1) The effective exhaust velocity ¢ . Unless otherwise 
indiccted all values of the effective exheust velocity referred 


to in this analysis are those definec by the equation 
C= constant /Ah 


where Mh is the enthalpy change for an isentropic expinsion., 
Because of the agsuaptions made these ere theoretical meximun 
values. In practical design work, en experimental velue of ¢ 

is obtained by dividing the mezsured value of thrust by the mass 
flow-rate, This experimental value includes pressure, friction, 
and divergence effects, end is sonevhat less than the value de- 


fined above. 


(2) The specific impulse Isp The specific impulse of a 
propellent muy be defined as the thrust per unit weight-rate of 
flow or the pounds of thrust per pound of propellent burned per 
seconc. It is obtainec from the equation | 


Isp = e/g 


. 


where g is the acceleration of gvevity, 
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PART II; wAMPLF CALCULATION FORK AN OVER-OXIDIZED uIXTURE 
OF LIGUID AMMONIA AND LIGUID CZCRE | | 


The objective of this section is to present in a logical 
sequence the c&leuletion of the data required to construct an 
enthalpy-entropy diagram and the calculations involved in 
applying the diagrams to the determination of rocket motor per- 
formance parameters. 

The products of reaction resulting from the combustion of 
LOS grams of liquid ozone and 5:,.0945 grains of liguié amzonia 
are considered. This is an over-oxidized system ( mixture ratio 
1.v86 ), and it is chosen because the minor component values for 
this case are more sensitive to changes in the estimated quan- 
tities than for stoichiometric and under-oxidized systems. This 
sensitivity, and methods for solving the resulting complicat- 
ions, are discussed, 

Data for the point on tne h-s diagram corresrcnding to a 
tempereture of ofC00O°K and a pressure of 500 PSIA ( 26.408 
atmospheres ) is obtained.’ In compiling the thermodynamic data 
from which to construct the charts, minor component caleulotions 
hed already been made for this pressure at 28c0°"K, Z000°K, and 
S500°K, and from a plot cf the sensible heat of the products of 
reaction and the heat ivalluble from the reacticn versus temper- 
ature the udiabatic flame tempereture T, had been estimated to 


lie between 3f00°K and SdC0°C*K. 


ofeP A: WINGH COMPONANT CALCULATICNS, 


The object of this step is to obtain the ecuilibrium com- 
position of the products of resection at S£Z00°K end 300 PSIa 
with the value of each component correct to the nearest theus- 
anth of a mole, The e:uations used are those developed by Dr. 
David Altman and Dr. Sidney Weinbaum of the Jet Propulsion 
Laboratory. This development is presented in Appendix I, 


The eyuation for this reecticn at elevated temperéetures is 


2.9412 WH, + 2.0835 O4 we Ho + bicO + eOg + PN + gNC + HOH 


+i4H + jJ0 


where the reactents are in the liguid state, 


From Appendix I the equations te be solved are: 


ip) baw B- a+ Kya’ 





oa oy 
at 
(ri) gm) Otte | 1, 3g. éKeb 
aé e 
cust 
(IIT) fi i 4+ g-H ow F(III) = 0 
Ee. 
D- 
2' | - 
c= be % f = 3(N = g) ie BX, 
a’ 2% 


J~2Ss/ 


ee 


Ate equilibriws constcnts in the ubove equations are’ those 
based on the ratios of mole fractions of the components in the 
reaction products, 

From the reaction equation, the reectants represent a conm- 
bined weight of 150.094 grams, and the values for the number of 


gram atoms of hydrogen, oxygen, and nitrogen are; 


H = &,823 
GO wm 6.5320 
N = £€,941l2 


Annother equation which is used as a check for "g" is 


g a 


. * % 
re % 





Steps in the caleuletion ere us follows; 


(1) A value of "n", the total number of moles of the reaetion 
products at 3£00°K and 300 PSIA, is assumed. Unless some data 
or previous results are available cn which to base an estimate, 
this value of "n" must be &@ guess, Hoxever, the calculaticns 

ere not too sensitive to errors in "n"®, and the second or third 
estimate based on the value obteinear in previcus trials should 
be close enough. The calculations are repeuted until the derived 
value egrees with the estimetea value to the desired degree of 
accuracy. In this case, component calculations fer 25cC0°K, 
S000°K, and 3500°K hed elreedy been made. A plot of "n" versus 


temperature wag prepered, and from this plot the value of "n" 


t f° 


at voLO°K wis estimetea te be 7.15, Using this value of mm", 
the velues of tiie necessary equilibrium evnstents based on the 
ratics of mole fractions gre determined from the tabuletea 
values of the corresponcing constants based on the ratios of 
partial pressures, 


The resulting values are: 


K-; .. 07859 


tt 


Ky m= , 000672 
Kao = 1605 
Kio sae 26454 


All numerieccl values in this calculation ere carried to at 
least four significent figures unless it is obvious tnat the de- 


sired zecuracy wili not be secrificed by carrying less, 


(2) A value of "a", the number of moles of Hy in the resetion 
products at 5500°K is assumed, This value must also be a guess 
unless previcus results ure available cn which to base un esti- 
mate, Here a plot of "a" versus temperature similar to that for 
"nN"! was prepared, and the value at 3200°K was estimated to be 


SOLv. 


(3) Using the above estimated values of "a" and "a", equations 
I, IL, and III are solved with the result that F(III) is equal 
to ~ £,1254, This negative value of F(III) indicates that .215 
4s too cmall « value of "a", From the variation of F(III) with 


Na" for previcus calculations, it is estimated that "a" must be 


a 


increased by about four thousunths, 


(4) The solution of equations I, II, and III is repeated for 
"ea" equal to . 19 and for "a" equal to .£20. The a 
values of F(III) obteined are - .<50'7 und .3886. Solving for 
all minor component Values for these two values of Nam, the 


results are; 


a rahe) eo 

b o. 8866 5, 6662 
a 01186 of TLS 
f 1.5949 1.4856 
ff baa. 17900 
h D545 . 5531 
of L779 0761 
j 1184 117g 


It is obvious from the cbove results that "g" and "e" are 
extremely sensitive to changes in "a", while the other compon- 
ents ere relatively ins emi Live, This extreme sensitivity is not 
typical cf stoichiometric ana under-oxidized systems, and when 
"a" nos been determined for these systems to the nearest thous- 
anth of a mole, other component values will also be correct to 
the seme cegree of ececurscy or se nearly so thet necessary cor- 
rections can be made by means of an stom balance check. For this 
czse, however, it is necessary to estimate a new value of "a" to 
four decimal places. This is uone by plotting F(III) versus "a" 
and lineerly interpclating for F(III) euusl to zero. A value of 


ie eel Go . oe LL nie cpt: . 


a 


(5) wocl¥ing @ weatiupe I, II, 4nd Ill *éth "2 osudi te yobea, 


l~z-rie ~ s wl > a ay ; ‘a « _ 7 = a8 ’ 
Fvl-i?} is equal to ony, e:nl other whlue’ ape Es Tolldtz,: 


& ee bot 
b 3, €v65 
e 07756 
g 1. ou 


a> & 


evoee 
i 0 ee 
- L162 


Cu. 


n 7267 


The derived value of "n" is now comparec with the estimated 
value. The diffcrence is .01¢<7. Since this is larger than the 
arbitrarily chosen maximum alicucble error cf .Cl, é new value 
of "n™ equal to 7.165 is G@zssumec, and the entire prccess is 


repested, 


(6) Using the velue 7.160 for "n", the eyuilibrium constants 


n 
Oo 


Up 


° 
6 


ie = Ly O7667 
3 
K = G8476 


= 166607 
= «6440 


~ 19 = 


(7) Since the error in the velue of "n" derived ebove was only 
slightly larger than thet aliowable, the component values will 
change very little, and it is still reasonably certain that the 
value of “at will lie between .£19 and sucOn In order Ge 
Gefinitely cetermine which of these ‘Ne the closer vélue, the 


calculations sre now repeated with "a" equal te .#195. 


(@) Using the value .£1285 for "a", the velue of F(III) is 

- ,O0.30. The value of "al" correct to the nearest thousanth of a 
mole is then .«:0. Using the check ec,uaticn fcr "g", the ecr- 
rect value is founc to be .160. Final results rounced to three 


decimal places ere: 


& <0 
b 3.680 
@ eS 
f£ Boece 
g LGD 
h 00% 
i 078 
d Lis 
n 7,166 


Bbexing an etom balance check: 


= ee 


‘ 


None of these values differs more then .001 from the value 
as determinea from the reactants. Also, "n! checks exactly with 
the assumed value, and all component values cen now be considered 


accurate to the desired degree, 


STEP Bs; CALCULATION OF PHE ADIABATIC FLARE TEEPLHATURE Tee 


when the procedure as outlinec in step A is carried out 


for S300°K, the following results ere obtainec: 


S200°K 2300°K 
@ e220 286 
bd 2.486 6.751 
e 0776 767 
£ 1.091 1.585 
gE 160 172 
h se4 5d 
i .078 116 
j wt ~1L59 
n 7,166 Toe 


The adiabatic flame teupercture qT. is defined as that teui- 
perature at which the sensible heet cf the ecubustion reaction 
products is equal te the heet released by the reaction at so0°K, 


The sensible heat is vefinec by: 


' ied Ci mC. 
y Tva >*- 
& C7 ; por — » z 
Sensible Heat = a Diy BO0¢K 


where "nj" is the number of moles of the i,, component of tia 


ya 
*.s 
I on i 


combusticn reaction products at the tempercture indicetea by 
the superscript ( “i we HEN, ee, FE" Bc. )» and "Ou," is 
the sensible hezt of the din component of the reacticn products 
at the temperature indicated by the superscript and as referred 
to the defined reference state. The values of "OH," as defined 
here are the enthelpy values tabulated in Table II. 


The heat released by the combustion reection is defined 


by: 
POR inl 
fev = : m Or. - YY (Reactants) 


where "ery" is the heat of formation of the 41,, compenent of 


the reaction products, 
In this case: 
g200°R . az00°K 
ny Hisoo°k = 220K 22,927 + 3,866 % 55,006 
+.776 X25,344 +1.391 25,947 


+.160 X24,519 + 534 X25,227 
+ 078 X14.405 + 116 14,405 


~05,45¢ keals per 7,165 moles 


il 


or per 150.084 grams. 


3300°K 3300°K a ae am ' . 
> ni Hy ayo = 1886 X 25,618 + 3,751 %54,362 


{ 


+ .767 *£6.308 + 1,585 24,837 

+ 17 %2L5.421 * .654%*% 24,121 

+ 116 *14,.901 + .159%*14,901 
£14,044 keals per 7.£70 moles 


or per 150,094 grams. 


me 
om fhe 


eo anOry 
on OO K 
t 


5,866 X 57.798 + ,160 *(-21.6) + ,554 *(-10.06) 


“ay 
+ .078 x (-52.083) + .118 * (-59,156¢) 
~ £.9412% 16,07 - £,0884 X (-31,6) 
= «23.090 keals per 7.165 moles 
or per 150.094 grams, 
S500°K —e 
he = 5.751% 57.798 +.172 « (-21,6) + .634 *(-10.06) 


+ .116 * (-52.089) + .159 x (-59.153) 


= £09.617 keals per 7.270 moles 


or per 150,004 grams. 


Using a linear interpolaticn to find at what temperature 
the heet released is equal to the sensible hest of the reaction 


products, the result is 


T = 3280°K. 


$ 
ey) 
CX 
j 


STEP C: CALCULATION OF ENTHALPY VaLuss, 


The enthalpy of a gus mixture at a temperature of T°% 
whose reference state is the ejuilibrium ccuposition at 400°K 


Ms, 


is given by: 
eo 7 : rae T°K 2 SOOTK _ at 
LB Hes, . OOK = Oy DEG s00 Oy ae (a; 3 Rey 


where all symbols wre the sume as previously cefined, 


3200°K 


Han 9% 205,432 +(4.412 ~ 5,886) * 57.798 


~ .160 * (-21,6) ~ 554 *(-10.06) 
- .O76 ¥(-52,083) - .116 *(-59.159) 


= 255.705 keels per 7.165 moles 
or per 150.084 grams. 


The enthalpy value plotted is cefined by: 


T°K TOK 
Dbz 9% = AHso0% / ® 


where "m" is the number cf grams of propellent mixture, in 


this case 159.004 grams. 


TOK 5200°K 
| ey cA | - 856.705 
5 Bs 00° nis Me 0°x 150,094 


1.704 keals per gran. 


it 





Pressure lines on the h-s diagram whose enthalpy values 
are calculated using Die oe as defined above correspond to 
equilibrium flow because the value hg ccntinucusly chenges 
as the equilibrium chifts. If pressure lings corresponding to 


constant composition flow are desired, it is necessary only to 


‘ue 4 -* TOK 
replace n, by n; in the equation for DH. 50% . The 
second summation term beconcs a constant value which is always 


Se 
greater than the value when ny * is used, The enthalpy value 


for a4 temperature of T°K end ccrresponding tc constant compos- 
jtion flow is therefore larger than the velue at T°K corres- 
ponding to equilibrium flow. Because the heats of fcormetion of 
the minor cemponents are negative, the fact that the minor con- 
ponents do not disappear in a constant composition flow means 
that less of the totul heat present in the gas mixture during 
such a flow is sensible, sit os a result T, and the exhaust 


velocity are both less than the velues for equilibrium flow, 
It should elso be noted that when equilibrium flow is 
being considered and enthalpy values are being calculated for 
temperatures less than 2000°K, the second summation term in 
the equation for Ate toe is zero because for these temper- 
atures the composition of the gas mixture is assumed to be the 


same as the equilibrium composition at S00°K, 


o~ ES 
oO —— 


oloP D: CALCULATION CF ENTROPY VALUES. 


The totel absolute entropy cf a ges mixture at 7T°K and an 


absolute pres»ure of "p% atmospheres is given by: 


: Tox G ed | 
e) e 
S ROK = ny Sy TK - on Ro in p 


Gx, O 
«(Cw 2K 
~ > n, R, in (n;/0) 


where Sincx is the absolute entropy of the 1h component of 
the reaction products, ana R, is the universal gas constant. 
Values of the abdsolute entrcpy for the gas components are 
tebulated in Table III. For this report, the values of n,/n 
were carried to four decimal places, and the natural logs of 
these fractions were taken from Volume III of the tables con- 
piled by the Works Progress Administraticn, which lists the 
entering argument from .U001 to 5.0.00 in steps of ,0001. 


For 3£00°K and 20.408 atmospheres pressure} 


_ = 220% 49,04 + 5,6u6 * 69,01 + .776 X 68,57 
ea 


+ 1.591 *64,31 + .160*69.47 + .534%*61.97 
+ .078* 39.15 + .118%* 49.79 


~ 7.163% 1,986 in 26,408 ~ 1,986] .226 In .220/7.16% 
+ 3,866 In 3,686/7.163 + .776 ln .776/7.16% 

+1,591 In 1.481/7.160 + ,160 In .160/7.163 

+ 534 In ,564/7.165 + .078 ln .078/7.163 

+ 116 ln .118/7.163 


- £6 


474.762 - 42,898 L¥.820 


D 
3200°K 
= 451.604 cals per degree Kelvin per 7.163 


moles or per 150,094 grams, 


The entropy value plotted is given by; 


For 3£00°K and 300 PSIA:. 


© 
so... = 45).634 
o200°K 150,094 


6.009 cals per gram per degree Kelvin. 


To calculate entropy values for pressure lines correspon- 


ding to constant composition flow it is necessary only to re- 
e) 


Tk Te fo) 
place ny by ny; in the equation for 8 0x , 


STEP Es CALCULATICN OF EFFECTIVE EAHAUST VeELOCITY and 
SPECIFIC IMPULSE, 


The effective exhaust velocity is given by: 





cm per second 


where "J" is the mechanical ecuivalent of heat equal to 


4,186 se i ergs per kcal, 


cg 


Since the plotted enthalpy value is OH/m the equation 


for "c" can be written as follows: 
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= 
cx g dn % 4,186 * 10 cm per second 


© 


The numerical faetor .0S261 converts cm per second to 


feet per second, The equation for "c" then becomes: 


2 eo 2 3 ae ie 
¢ m= ,06281 X 2.69545 x 10 ak, fect per second 
€ 


Finally: 


/ if 
ec = £9435.41 An, _ Leet per second 
) 


m 
The value of 4h,° is very simply determined from the 
2 


h-s diagram. A verticél line is dropped fron qT. on the pressure 
line corresponding to the chamber pressure to the pressure line 


corresponding to the ambient exhaust pressure, The value of en- 


€ 


vote ° The enthalpy 


z 
thelpy for this latter intercept is Odin 


crop is then given by: 


i 


i T Z 
An® = Ah - At®.. 
ae 500°K 300°K 


For the sumple calculation the isentropic expansion fron 


400 PSIA to 14.7 PSIA is considered, For this cése: 


=i 


a5 . , 
Oho — 1.8:1 - 1.002 m= .799 keals per gram 


Le 


from Figure <2, 


0493.41. / .79Y feet per second 


ex 
c = 6486 feet per second’ 
Isp is vetermined by cividing the value of the effective 


exhaust velocity by the value of the acceleration of gravity. 
For this culculaticn: 


Ten = — — <64 seconds 


STEP PF; DETERMINATION OF THs EXHAUST TEMP:RETURL . ° 


The value of the exhaust tempereture qT. can be cetermined 
roughly from the h-s diagram. It is the temperature value at 
the intercept of the isentropic expansion line and the pressure 
line representing the ambient exhaust pressure. A more accurate 
value is obtained from a plot of entropy versus temperature, 
and all values reported in this anelysis are determined in this 
manner, Cn the T-s diagram a vertical line ig dropped from T, 
on the pressure line representing the chamber pressure to the 
pressure line corresponding to the ambient exhaust pressure. 


fae] 


The temperature value of this intercept is os 


PART III: DISCUSSION OF RESULTS. 


There is no basis for «n extended discussion of the per- 
formence values reported in this analysis. No experimental 
data for the propellent systems considered here was available, 
and the only analytically determined values available for com- 
parison were those for the stoichiometric mixture of liquid 
ammonia and HFNA calculated by the Jet Propulsion Léboratory. 
The results obtained by — methods sre as follows: 


A, Equilibrium Flow; 


arameter ri Valu Gra al Value 
f. 2600°K £599°K 
- 1631°K 1640°K 
c 7815 ft/sec 7324 ft/sec 
Isp PETE see 227, £608 


¥, 


B. Constant Composition Flow: 


Parameter JPL Value Graphical Value 


wy S600°K 2590°K 

2 1505°K 1508°K 

c T17 ft/sec 71isy ft/sec 
Le co&.l sec £25,8 sec 


It should be noted that ts is determined in the same way 
for both analytical and graphieal solutions. For the graphical 
solution, however, T, is reac directly from the eppropriate 


T-s dlagram, and the values of "ce" and I depend only on the 


5 
value of the enthalpy chenge for the en considered, The 
vélue of this enthalpy change is tuken directly from the h-s 
aLagran, The fact that the above values egree to within two 
tenths of one percent is considered 2 positive check on the 
accuracy of the graphical methed. 

Pressure lines on the h-s and T-s diagrams corresponding 
tc constant composition flow are constructed cnly for the ex- 
pansicn from a chamber pressure of 5090 PSIA to an ambient’ 
exhaust pressure of one atmosphere, These pressure lines are 
indicated on the dlegrans by dashed lines, and should be used 
only for the above specified expansion. The error, however, 
Will not be great if they are used to epproximate the con- 
stant composition flew paranecter values for 4& chamber press- 
ure of 60C PSIA, The fact that exhaust pressure lines for 
constant compesition flow are included for only this one con- 
dition is not intended to convey the impression thet the 
actual expansions of these propellent systems conform to the 
assumptions for equilibrium flow. It was originally: intended 
to calculate complete constant compesition flow data, but. 
time limitations prevented this. In the absence cf more com- 
plete date, no definite conclusions can be dréwn as to which 


type of flow would result from the combustion of the propell- 


ent systems considered, 


For very hot systems, such as HC, » the position of the 
exhaust pressure lines representing constant composition flow 
on the h-s and T-s diugrins is cuite different from that of the 
ccrresponcing pressure lines for equilibrium flow. This dis-— 
plecement on a percentage busis is greatest on the T-s diagram, 
and <s a result T Ls the performance parameter whose equili- 
brium flov value differs most from its constant composition 
flow value. For ths W.-0, system the T, value for constant 
compositicn flow is several hunured degrees lower then the 
equilibrium ficw value, 

The plots of specific impulse versus mixture ratio show 
one significant Sau eaey. Toe mixture riutic for which the 
maximum value of the specific iwpulse is obtained appsarentiy 
shifts towards the over-oxidized side us the exhaust pressure 
approaches absolute zero. If this effeet is general, it will 


be an important factor in the design of the latter stages of 


multiple-step rockets. 


#2 ; 


PART IV: CONCLUSIONS AND RECOKKLNDATIONS, 


Because of the extensive calculations which are required 
to construct even a sxeleton enthalpy-entropy. diagrai for a 
propellent system, this graphical wethod is not recommended 
unless & comprehensive study of the system is to be made, It 
does, however, have several cefinite advantages over analyti- 
cal methods. The h-s diagram makes possible the presentation 
of a large amount of data in a coneise und readily useable 
form. If the dlagram is a fairly complete one, an almost infi- 
nite variaticn of the physical boundary conditicns defining 
the expansion process is possible, whereas one set of analyti- 
cal calculations is limited to one combusticn concition and 
one exhaust condition, | 
| The graphical method does not require that attention be 
focused on either the assumptions of equilibrium flow or of 
constént composition flow to the exclusion of the other since 
exhaust pressure lines for both conditions can be plotted on 
the same diagram. Neither is it nétes sey to assume isentropic 
flow if some criteria for estimating a "condition curve” are 
aveilable. 

In order to better evaluate the usefulness of the h-s 
and T-s diagrams in cCilculating perfcrmance perameters, this 
graphical method should be applied te & propellent dystou tor 
which complete experimental data is available. If @ complete 


set of pressure and temperiture lines are constructed, not 


=e” « 


only for the thermodynemic conditicns in the rocket motcr 
chamber and at the nozzle exhuust, but else for conditions 
near the throat, the h-s cisgrzm could be applied to the cal- — 
culetion of the characteristic velocity and the nozzle thrust 
coefficient as well as the parameters determined in this 

analysis. It is believed that a therough investigation of one 


propellent system in this manner would be of greét value, 
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APPENDIZA I 


MINOR COMPONLRNT CALCULATIGNS OF 


6, H, 0, N SYsTexs 


The most tedious part cf a performance eczlevuletion for 


high temperature systems is the solution for the gas composition 


at chamber and exit conditions. For systems containing carbon, 


hycrogen, oxygen, and nitrogen, there are ten components to 


sclve for if the chamber temperature exceeds about 2500°K and 


an eleventh (atomic nitrogen) if. the temperature exceeds about 


3000°K, Trial and error solutions require many cycles of cper- 


ations and so are fairly long. it is possible, however, to set 


up ten equaticns involving the ten’‘unknown components in terns 


of the equilibrium constants and the mass balance equations. 


These equations can then be reduced in number by eliminating 


unknowns end the resulting equations cen be solved directly 


for the components. One such scheme is the following: 


Let the number of moles of the components be represented 


py the following symbols; 


a 


b 


These ten unkmowns are 


—— 
Amal 


i = % : 
Be, 
6 = yo 
h = Noy 
i= Dey 
= Nn 
J © 


reluted by the fcllowing ecuetions: 


a 


@ B= (6) hac = Kyo 
(©) ga = K, (7) e-d=- Cc 
(3) ie * Ke (8) O-c¢- Saag se AVe Go 
(4) ai = hey (9) 4-12 .% = H 
(5) i= K, (io) Sf-g = & 
a® 


where C, 0, H, ona N represont the total number of gram etoms 
of the element in the sample chosen for the czlculation,. The 
K's are given in terms of the number of moles of the consti- 
tuents and are related to the K's (in terms of pressures) by 


the relation: 
nA. - 
K = K,(n/p) P R 


where n is the total number of moles, P, the total pressure, 
and Ny-Np represents the difference in the number of moles 
between the products and reactants for the particular reaction 
under consideration, | 

The elimination of the unknowms may now be accomplished 
«aS Yoliows. Solve Eqn. for c:3 

C= C- a 
Substitute into Econ. (1) to yield 
| bt - bd - Kjad = 0 


This gives for a: 


d= bC 
b-K 





Le 


Substitute back into Eqn. (7) and obtein for e: 


e .~. C-_ bf = Cc 
b- Kya i - ae 
ha 


The remaining unxAnowns wre now readily solved for in terms of 


a@ and b to yield: 


b - Aya \ as 
& *. — 
e - _b Ke iif: aa ‘g 
_” . 
f =- $(it-g) 4 = _» Ke 
2. 


Substitute for h end Z into Eqn. (3) and obtain for pb: 


H- 2a - Koa" 





(I) b = : 
 & R10 


ai 





Solve for g from Eqn. (8) to yield: 











Kia — xb h10 K 
@ piano .c¢ a iSene tone 
Ke - b ai a 

a ae 

= 2h, b 


Another equation for g is obtained by elininating f from 


Eqn's. (£) and (10) to give: 
@ 
(IIT) hm - g 2h = © @ P(TIT) 


The equations I, Ii, and Iii form the besis for the scluticn 


of the unknowns. This is accomplished us follows: » 


PROCEDURE 
(1) Choese a likely value of a, 
(2). Caleulate the corresponding value of b from (I). 
(4) Caleulate the corresponding value of g from (II). 
(4) Substitute the values of a, b, and g into qn. (III) 
to obtain F(III). 


(5) Hepeat ealculetions and interpolate to those vélues 


which make F(III) s usl to zero. 


in general, it will be found that if F(III) is greater 


then zero, & should be decreased, and vice versa, 


The following are special céses: 


A, Case where there is no C in the system 2%, Cog 


Eqn's, (I) and (III) are the same but (II) beeomes: 





t= c8fo2-o --— > 
= 
cB ak c 


€3 


The procedure is the sane. 


B. Case where there is no N in the system C, H, O ; 


Eon. (I) is the seme, (III) is elimineted, and (II) becomes; 











K K oy K,a =— 2b 
io ”" . EKGD we ie 
(II) b od - -  —————— 
a2 . a Ky,a - b 


The proper & will make F(II) eyual to zero. 


C, Case where both C and N ure absent H, © ; 


Eon. (I) 4s the same, (III) is eliminated, and (II) becomes: 


—sk Ge = 





K K migb 
a’ " & 


The proper & makes F(II) equal to zero. 
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FIGURE 18 


THEORETICAL SPECIFIC IMPULSE FROM h-s DIAGRAM 
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TASLE IV 


HEATS OF FORMATION AT CONSTANT PRESSURE 
USED IN THIS ANALYSIS 


Heat of Formation in 


Constituent Kilocalories per Mole 
H»0 +57.798 
NO ~21.600 
OH -~10.060 
H -52.089 
0 59-159 
NH ( 1) +16.070 
NOo ~ 6.100 
HNOs +41 .660 
03(1) -31.500 
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ATOMIC WEIGHTS USED 


Blement Atomic Weight 
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1.008 
16.000 
Ho 2.016 
Ho0 34.016 
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(DD ble 00? 2.082 £6882 6.090) 

ae 153i ‘1100 161] 73 : 
| .007037, 0234), .0289 


02595, 025% 
“ios owsheg 1 socenae — 
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TABLE YL 
SAMPLE CALCULATION SHEET (A) 


P=300 PSIA 
T=3200°K 

















ee er BALL 
| jyoloeaa4l | bh} 
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ave | 4626] .4690| .4690| 4634|.4685| iN 12.992 
J 0462| 0476) .0484| ofl ose 
2. 2 
Bihe I37Z. 1213730213746 
ie "0710 079.0781 es 
B.3166) 3.307718, 3055 |8,3068'9.3065| | 
ae Sa AEE ere 
be NSE 15,1057 (81026 [91047 151018 _ 
ZKab/at | A1563, 4007 24702991 99G 
03103 .0305 ines 0304 
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TABLE WIL 
SAMPLE CALCULATION SHEET (B) 


NH, — Oz (m = 150.094 GRAMS) P= 300 PSIA 


OVER-OXIDIZED (r =1.996) Q. (REACTANT S)=—I8. 3589 KCALS 


ES SS A i NS Sates 


eT Tose [ie [ane [ee [see] ese] — 

—— 
aie | ata ae aaa ae a 
‘iA OR 











4a | _11.438 (.400//.301 1. 386 11,385 1 3oe 
eee eee ao 

nou | toa ato! 534 | aie | ese 53 
| 032 242 


0 |_| ‘008 | 063 | tia [sy | e592 
ee ee | | | ee ee 
| on | | @ 84 | 7,009 | 7.103 | 7.249 | 7.270 17.554 
_ “ene mas?) 
IZnjAWi | __—(/48,/401/88,584|z05.432, 214,044 231,110 
| | | «deel 
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TABLE VIIT 


ENTHALPY - ENTROPY SUMMARY 


P = pressure in psia 


Oh = keal per gran 


s = cal per gram per °K 





RFNA - NHz OVER-OXIDIZED (r = 3.00) 























Temp°K P=» 600 Ps 300 Pwilkt.7 Pw 7.35 =P» .147 
3000 Ah 1.7381 1.4202 
+s | 2.811 2.886 
2500 Ah 1.004 1.012 
s 2.674 2.739 
2400 Ah 0.943 0.948 
s 2.649 2.713 
2300 Ah 0.885 0.888 
$s 2.624 2.688 
2000 Ah 0.723 0.724 0.730 
s 2.549 2.611 2.884 
1500 Ah 0.479 0.479 0.479 
8 2.740 2.802 3.152 
1300 Ah 0.390 0.390 
: 8 2.677 2.739 
1000 Ah 0.261 0.261 0.261 
8 2.564 2.626 2.976 
500 ah 0.070 0.070 
8 2.362 2.712 
a0 Ah 0.035 
8 2.634 
300 Sh 0.0 0.0 
s 2.185 e535 


ma 8 6S% 
eB ‘wee 
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RFNA - WH, STOICHIOMETRIC ( r = 2.21) 














Temp K P= 600 P= 300 P#ilh.7 P#7-35 P #147 
3000 Ah 1.517 1.57 
8 ee 946 53.0353 
2T00 Ah 1.233 
s 2.846 
2600 Ah 1.154 1.171 
































2000 Ah 06.770 0.772 0.783 
a 2.650 2.716 3.005 
1200 Ah - 0.610 
9 2.911 
1500; Ah 0.510 0.510 0.510 
5 2.849 2.914 3.282 
1400 Ah 0.462 
5 2.815 
1300 Ah 0.414 
rm 2.885 
1000 Ah 0.278 0.278 0.277 
8 2.661 2.726 3.094 
500 Ah 0.074% 
8 2.813 
300 Ah 


~~ 


RFNA - NHs UND@R-OXIDIZED ( r= 1.50 ) 





femp°K P= 600 P= 300 P#14.7 P= 7.35 P= 147 





3000 Ah 1.548 1.583 
s 3-272 3.362 





2500 ah 1.169 1-172 
es 3.134 3.213 





2100 Ah 0-920 0.921 
s 3.026 3.103 





2000 Ah 0.8641 0.861 0.864 
s 2.997 3.074 3.413 




















1500 4h 0.577 0.577 0.577 
8 Be 248 Je 326 3.764 
1400 ah 0.418 
8 3.130 
1000 4h 0.317 0.317 0.317 
$ 3.038 3.116 56655 
900 Ah 0.268 0.268 
800 Ah 0.221 
8 3.008 
500 Ah | 0.086 0.086 
a 2.797 3.234 
400 Ah 0.042 
8 3.138 
300 ah 


@ 
\W O 
@ e 
Oo 
a) 
~“ 





3500 Ah 2.042 2.208 
s 3.034 %.160 





3400 Ah 1.888 
gs 2.999 


3300 Ah 1.749 1.850 
S 2.958 De 655 





3200 Anh 1.704 
8 3.009 





3000 46h 1.415 1.458 1.940 
s 2.852 2.930 3.379 

















2500 Ah 1.034 1.107 
s 2.776 3.080 
2000 ah 0.747 0.747 
8 2.920 2.984 
1590 ah 0.488 0.488 0.488 
Ss 2.773 2.835 3.187 
19090 4h 0.266 
8 3.008 
500 Ab ° 0.071 


3 2.739 


on T 13 = 
NHz - O5 STOICHIOMETRIC (r = 1.41) 


Temp °K P= 600 Pe 300 Peil4.7 P=7-35 P= .i47 


a. 





3500 Ah 2.351 2.554 
s 3.278 3.412 


3400 Ah 2.167 
s 3.225 


33900 Ah 2.9001 2.) 
$ 3.175 3.288 

















3200) 4h 1.956 
g 3.233 
3000 Ah 1.597 1-658 2.267 . 
Ss De 046 32137 3.657 
2500 Ah 1.146 1.260 le 5e) 
8 2.951] De 2904 3 383 
2000 An 9.819 0.824 
Fs) 3+ 099 3.169 
1800 Ah 0.692. 0.693 
3 3-031 3.098 
1500 4h 0.555 0.533 0.533 
3 2.934 7.001 3. 380 
1000 ah 0.289 
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WH, - 0, UNDER-OXIDIZED (r = 1.10) 



































Temp°K P= 600 P= 300 Pp #14.7 P=7.35 P= .147 
3500 Ah 2.275 2.469 
ag 3.508 3.647 
3200 Ah 1.810 1.892 
s 3.369 3.474 
3100 Ah 1.692 1.750 
s 3.332 3.429 
3000 Ah 1.588 1.628 2.164 
s 3.297 3. 389 3-919 
2500 Ah 1.195 1.249 
S 3.232 3.590 
2000 Ah 0.878 0.879 
a 3.426 3.504 
1700 Ah 0.699 0.699 
2 3.329 3.406 
1500 Ah 0.585 0.585 0.585 
8 3.858 3-335 3.770 
1000 Ah 0.320 
s 3.556 
500 Ab 0.086 


3-233 


=" a2 - 


NHs ~ Oz UND@R-OXIDIZED (r = 0.80) 

















Temp°K P= 600 P= 300 P#14.7 P= 7.35 P= 147 
3500 Ah 2.563 
3 3.994 
3000 Ah 1.754 1.772 
8 3.660 3-752 
2600 Ah 391 1.398 


le 
s 3.526 3.619 


2500 4h 1.317 1.322 1.377 
8 3.497 De 589 4, 004% 





2000 Ah 0.975 
8 3.826 
1500 Ah 0.653 0.653 0.653 
1400 Ah 0.592 
3 3-599 
1300 Ah 0.532 0.533 
8 3-555 3.648 
1200 Ah 0.474 0.474 
& 3508 3-599 
1100 ah 0.417 
8 3.548 
1000 4h 0. 360 0.360 0.360 
a 3.404 3.494 4.002 
500 Ah 0.098 


8. 3.640 
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H, - Os OVER-OXIDIZED (r = 9.92) 





Temp°?xK P2600 P# 300 Pwlk.7 Pe 7.35 P= .147 
8OOO) = 6 Ah 3.935 #555 
s 3-839 & .087 





3700 Ah 2.993 3-350 
s 3.595 3.774 


3600 Aan 3.028 
g 3.686 


3500 Ah 2.516 2.746 5 B45 
8 5 462 36 607 4 © T85 











3000 Ah 1.706 1.765 2-431 2.777 
s 3.213 3.306 3.858 4.055 
2500 ah 1.334 1+ 383 2.367 
3.464 L Phe 4. 380 
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Hp - Os STOICHIOMETRIC (r = 7.94) 








Temp°K P «600 P= 300 Pewlk.7 Pe« 7.35 Ps .1A7 
4000 Ah 4.580 5.293 
s 4.217 49 
3800 Ah 3.819 
s 4.023 
3700 ah 3.495 3.912 
s %3.936 A.1AL 
3600 ah 3.542 
8 4.037 
3500 Ah 2.936 3.213 6.204 
+. 3.780 3.947 5.246 
3000 Oh 1.965 2.046 2.869 3.278 
= 3.581 3.588 4.229 b 457 
2500 Ah 1.397 1.589 1.617 2.825 
8 3.352 3.752 3.859 4.821 
2000 Ah 0.992 0.999 1.092 
| 3.507 3.587 4.069 
1800 Ah 0.833 
4 3.930 
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No - 03 UNDER-OXIDI ZED (vr = 6.61) 





Temp°?K P= 600 P# 300 Pwril4.7 P=#T7.35 P= .1L47 








4000 Ah 4.680 5.454 
s 4.600 4.910 





3800 ah 3.861 




















s 4.392 
3700 Qn 3.516 3.957 
s 4.298 4.522 
3600 Ah 3.561 
Ss 4A] 
3500 Ahn 2.933 3.213 6.460 
s 4.136 4.314 5.738 
3000 Ah 2.000 2.058 2.836 3.267 
s 3.850 3.960 4.632 4.880 
2500 ah 1.490 1.565 1.611 2.778 
5 3.754 @8177 4.286 5.294 
2000 Ah 1.086 1.090 1.124 
5 3.966 4.057 4.582 
2500 Ah 9.719 
8 &. 350 
1000 Ah 0. 391 
8 4.084 


=~ & 19 ~ 
Ho - Os, UNDER-OXIDIZAD (r = 4.96) 


Temp °K PeG00 P= 300 Peilt.7 Pe 7.35 P= .147 


aes 





4900 Ah 5.108 5.949 
s 5.367 D7 1 TF 











3600 ah 3.569 
s 4.960 
3500 Ah 3.299 %.560 T-OT2 
be 4.886 5.083 6.690 
2400 8 8=6ah 3.256 
g 4.994 
4000 Ah 2.359 2.41) 3.125 3.556 
“gy 4.598 4.732 5.490 5.760 
2500 Oh 1.781 1.863 1.908 3.023 
8 4.503 SO, 5.171 6. 297 
2000 oh 1.3597 1.408 i es 
S 4.792 4.907 5.579 
1500 An 0.870 
8 5. 301 
1000 Anh 0.77. 
3 4.983 
800 Ah 0. 332 


= ae 
TABL™ IX 


PERFORMANCE SUMMARY 
RFWA - NH3(1) 


ee ee ee a 
r Pe Pe C Isp To K Te K 
3.00 300 14.7 6645 207 2325 1360 
2.21 (stoic.) 300 14.7 7323 228 2599 1630 
1.50 300 14.7 6898 214 2062 1110 
3.00 600 14.7 Tl 223 2331 1170 
2.21 (stoic.) 600 14.7 79K2 uy 2621 1425 
1.50 600 14.7 7384 230 2063 955 





3.00 300 7.35 7167 225 2325 1190 
2.21 (stoic.) 300 7.35 7914 246 2599 1425 

»50 300 7.35 7384 230 2062 955 
%.00 600 7.35 7576 236 2331 1010 
2.21 (stoic) 600 7.35 8437 262 2621 1240 
1.50 600 Tao5 7782 eho 2063 810 
3.00 300 0147 8721 271 2325 480 


2.21 (stoic) 300 147 9751 303 2599 615 
1.50 300 »147 8869 276 2062 350 


600 147 8854 275 2331 400 
(stoic.) 600 0147 9920 308 2622 515 
600 -147 8996 280 2063 300 


ee 
oNs 
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TABLE IX 


PERFORMANCE SUMMARY 


a eS ee 
r Po Po C Zap To K Te K 
2.00 300 14.7 84B6 264 3280 2400 
1.41 (stoic.) 300 14.7 8930 278 3285 2460 
1.10 300 14.7 8902 277 3124 2050 
0.80 300 14.7 8292 258 2564 1423 
2.00 600 14.7 9252 287 3352 2200 
1.41 (stoto) 600 14.7 9749 303 3362 2290 
1.10 600 14.7 9649 300 3179 1802 
0.80 600 14.7 8910 277 2573 1220 
2.00 300 7.35 9199 286 3280 2200 
1.41 (stoic.) 300 7.35 9692 301 3085 2280 
1.10 300 7.35 9621 299 2124 1805 
0.80 300 7.35 8921 277 2564 1230 
2.00 600 7.35 987 307 3352 2185 
1.42 (stoic.) 600 7.35 10414 324 3362 2100 
1.10 600 7.35 10233 318 3179 1580 
0.80 600 7.35 9hl2 293 2573 1050 
2.00 300 -147 =611 684 363 3280 1100 
1.41 (stoic.) 300 o147 8 12409 386 3285 1220 
1.10 300 -147 #811872 369 3124 796 
0.80 300 147 10786 335 2564 465 
2.00 600 -147 11936 371 3352 1140 
1.41 (stoic.) 600 147 = =©12722 396 3362 1060 
1.10 600 e147 =. 12130 377 3179 
0.80 600 0147 = 0972 34) 2573 380 


- T 22 = 
TAME TZ 


PSRPORMANCS SUMMARY 


H(1) - 03(2) 


—— $$$ $$ 
Yr Po Pe Cc I 8D Te K Te o% 
9.92 300 14.7 10129 315 3634 2870 
7.94 (stoic.) 300 14.7 10613 330 3640 2900 
6.61 300 14.7 11107 8345 3607 2820 
4.96 300 14.7 11762 366 347 2500 
9.92 600 14.7 11059 344 3752 2770 
7-94 (stoic.) 600 14.7 11665 363 3760 2815 
6.61 600 14.7 12227 380 3721 2710 
4.96 600 14.7 12816 398 3537 2280 
9.92 300 735 11021 343 3634 2735 
7.94 (stoic.) 4300 7.35 11629 362 3640 2765 
6.61 300 7.35 WAP 877 3607 2670 
4.96 300 7235 12768 397 3447 2275 
9.92 600 7.35 11856 368 3752 2625 
7.94 (stoic.) 600 7.34 12521 389 3760 2675 
6.61 600 735 13081 407 3721 2550 
4.96 600 7.35 13651 424 3537 2010 





’ 300 147 18391 4&7 3634 1920 
94 (stoic.) 300 L147 15217 473 3640 2015 
300 aa 15768 490 3607 1650 
300 147 16081 500 34k7 1080 


600 0147 14819 461 3752 1710 
(stoic) 600 0147 15739 §=489 3760 1910 
600 147 16214 504 3721 1430 
600 .147 16442 511 3537 910 


rH A) 
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~ T 23 - 
TASLE X 


FROZEN PERFORMANCE SUMMARY 


RFNA - NH3 
r Po Pe C Isp To Ke K 


3.00 300 14.7 6563 204 2325 £1315 
2.21 (stoic.) 300 14.7 7199 224% 2599 1505 
1.50 300 14.7 6898 214% 2062 #£«®1110 


HH; - 05 


Be 300 14.7 8155 254% 3200 1882 
1.41 (stoic.) 300 14.7 8544 266 3285 1900 





1.10 300 14.7 8643 269 3124 1780 
0.80 300 14.7 9232 256 2564 1395 
Ho - Oz 


9.92 300 14.7 9578 298 3634 2095 
7.94% (stoic.) 300 14 10091 314 3640 #2100 
6.61 300 14.7 10566 329 3721 #$=:2060 
4.96 300 14.7 11306 352 3487 #1960 


~] 











fi2Se 


Thesis 
w555 White 
the calculation 


enthalpy-entroPy dia- 
grams for rocket pro~ 
pellant systems. 


of 


Sail 
; 


DUDLEY KNOX LIBRARY 


ha 





